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Abstract 

In this paper we suggest and theoretically study a tapered plasmonic nanostructure which con- 
nects the incident wave beam with a subwavelength spatial region where the field is locally enhanced 
in a broad frequency range or for different operation frequencies. This spatial region has a frequency 
stable location near the contour of the tapered structure. This results from a special waveguide 
mode which can also exist in the tapered structure. We foresee many possible applications for 
our structure from prospective near-field scanning optical microscopes to interconnects between 
conventional optical waveguides and prospective optical nanocircuits. 
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I. INTRODUCTION 



An actual problem for many branches of the modern optics is the concentration of the 
light energy initially transported by a conventional wave beam into a subwavelength spatial 
region. The problem is not new and there is a whole body of literature devoted to so-called 
subwavelength light concentrators which is difficult to refer in this paper. In most of these 
works the objective is to achieve very high resolution or very large local field enhancement at 
a certain frequency. In works on plasmonic nanostructures designed to enhance of thin-film 
solar cells (see, e.g., Ref. 1) the subwavelength field concentration is achieved together with 
efficient transmission of the light energy into this region (photovoltaic layer of nanometer 
thickness). However, as a rule, this subwavelength field concentration is a narrow-band effect. 
Rarely, it is achievable with the same nanostructure within 2-3 very narrow resonance bands 
(see, e.g., Ref. 2) 

Most sub-wavelength light concentrators such as a solid metal tip with nanoscale curva- 
ture of the apex or a nanoantenna operate in the regime when the surrounding space is also 
illuminated. On the contrary, tapered plasmonic nanostructures, we have unified under the 
name of metamaterial nanoptips (MMNT),-~- concentrate the light energy (in our simula- 
tions) in a needed sub-wavelength region so that the surrounding space is not illuminated. 
This way the parasitic excitation of the objects located outside the selected region with the 
incident wave beam with is avoided. Structures capable of realizing this regime are also not 
new in the nano-optics. They are actually being used in some schemes of field-enhanced 
nanosensing, especially of the field-enhanced fluorescence,- of the field-enhanced Raman 
spectroscopy- and also in field-enhanced microscopy 8 . Beside these applications, structures 
which realize such a regime could be used for interconnections between conventional optical 
information systems and prospective sub-wavelength devices. Sub-wavelength devices which 
need to be connected with usual optical waveguides are plasmonic nanowaveguides which 
have developed fast since the past decade,- 1 ^ optical nanocircuits,— and active all-optical 
nanodevices (nanoswitches, nanotransistors, nanodiodes and optical memory cells^). In all 
of these cases the general problem can be formulated as follows: One needs to transmit 
the optical signal from a conventional waveguide to a needed nano-objects whose locations 
are known while minimizing the cross-coupling to other nano-objects. This problem for- 
mulation recently led to the general idea of a specially shaped (e.g. tapered) waveguide in 
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which the usual wave beam formed by propagating spatial harmonics would concentrate in a 
sub- wavelength region outside the waveguide, for example, near the end of the waveguide.— 
We complement this problem formulation by a requirement of a broad frequency range or a 
multi-frequency operation. 

Notice that such a mode concentration was realized in some conventional tips used in 
scanning near-field optical microscopes (SNOM), namely in those being composed of a con- 
ically shaped optical fiber covered with a metal film having a nanometer aperture at the 
end. However, in such tips the hot spot is located inside the cone which causes almost all 
the incident wave to reflect from the end of the tip and the energy transmission through- 
out the tip is very low. Attempts to use the extraordinary transmission effectr^ 1 ^ for such 
structures are known. However, only a rather modest improvement of the transmittance was 
reported in these works and even this modest improvement was achieved by introduction 
of very complicated nanostructures, for example nanocorrugations in the metal film of the 
narrowest part of the nanocone.— 

More promising results were theoretically demonstrated for our MMNT. In Refs. 3 and 
17 tapered waveguides performed as metamaterial prisms (in Ref. 3) it comprised silver 
nanorods oriented along the prism edge, in Ref. 17 it was formed by silver nanoplates parallel 
to both the edge and the waveguide axis) efficiently transformed the incident plane wave 
illuminating the base of the wedge into a light nanojet. The light nanojet is a light beam 
keeping strongly sub- wavelength width (A/4) up to the distance A/2 from the nanotip apex. 
In Ref. 4 a pyramid filled with a metamaterial of plasmonic nanospheres (also illuminated 
from the base) allowed a very small hot spot (of the order A/100) at the apex. This light 
concentration was achieved in these works together with negligible scattering and acceptable 
reflection losses. 

Structures studied in the present paper possess important advantages compared to the 
pyramidal nanotip reported in Ref. 4. First, they are multi-resonant and the overlapping 
of resonances can result in a very broad operational frequency band. Secondly, the spatial 
domain of the enhanced electric field is nearly the same for all operation frequencies. The 
multi-frequency operation and the frequency stability advantageously distinguish our struc- 
ture from substantial plasmonic nanoparticles matched to a conventional optical waveguide 
in Ref. 19. We believe therefore that our MMNT reported will give a significant development 
to ideas of Ref. 19 and other relevant works published to the present time. 
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II. SUGGESTED STRUCTURE AND ITS OPERATION 



A. The main idea 

The proposed structure comprises parallel rectangular nanoparticles shaped as nanobars 
or as nanoplates of noble metal (silver or gold) separated from one another with dielectric 
spacers (also shaped as nanobars or nanoplates). Both variants are shown in Fig. Ufa) and 
(b). In the first case the structure width is much larger than its thickness. This MMNT can 
be located at the dielectric substrate and excited by an integrated (e.g. ribbon) dielectric 
waveguide having the same substrate. A rather similar multi-resonance structure has been 
known since 2008 and was called optical xylophone in Ref. 20. However, optical properties 
of the xylophone structure in which rather wide metal nanobars are separated with narrow 
spacers are different. The xylophone structure^ would better correspond to our Fig. Ufa) 
if we replace in this sketch metal and dielectric nanobars. The optical xylophone located 
at the dielectric interface was excited by an incident plane wave and at a given operation 
frequency the local field was enhanced near one metal nanoplate resonating at this frequency. 
Our structure (planar MMNT) is excited by the wave beam impinging the structure from 
the base of the structure, the field concentration holds in a frequency tangle spatila region, 
the resonances of different nanoplates overlap and form a very broad frequency range and 
all this becomes possible due to novel plasmonic waveguide modes that we called the edge 
modes. 

Another MMNT studied below is a pyramid performed of stacked metal nanoplates sep- 
arated with dielectric spacers Fig. [U(b). This nanostructure can be built at the end of an 
optical fiber. It seems to be similar to the plasmonic waveguide studied in Ref. 21. How- 
ever, due to different dimensions the operation of our structure is different from that of the 
waveguide structure in Ref. 21. The purpose of the study in Ref. 21 was the homogenization 
of the layered waveguide and the tapered variant of the waveguide was studied in Ref. 21 
only in order to demonstrate that the tapering is possible without destruction of the wave 
guidance. In Ref. 21 plasmon resonances of nanoplates were not exploited. In our study we 
combine these resonances together with the edge waveguide modes. 

In both planar and pyramidal sets of nanoparticles the size of the metal nanoparticles is 
gradually varied so that the plasmon resonances of the structure cover a sufficient part of 
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(a) (b) 

FIG. 1. (Color online.) Sketch of the proposed metamaterial structures for broadband local field 
enhancement with stable position of the hot spot: (a) - planar structure; (b) - pyramidal structure. 



the visible range. 

The idea of the surface and edge modes arose from the analysis of surface waves at the 
interface of the chopped medium with free space. We call as the chopped medium a periodic 
layered metal-dielectric structure with nanometer period. The sketch of this medium is 
shown in the left panel of Fig. |2j At the interface of the chopped medium with free space 
the surface wave excites which has a different dispersion than surface-plasmon polariton 
(SPP) supported by the interface of the solid metal. The dispersion of SPP for the solid 
silver sufficiently differs from the light line only within the frequency interval 780 — 970 THz. 
This interval refers to the ultraviolet range. In the visible frequency range the solid silver 
does not support surface waves. The SPP resonance of the chopped medium occurs in the 
visible range. In the example of the dispersion diagram presented in Fig. [2] silver layers are 
of thickness 5 nm and dielectric layers with e = 2.2 are of thickness 10 nm, whereas the 
SPP resonance holds at 780 THz. From 400 to 750 THz (the lower and upper bounds of the 
visible range) the slow-wave factor is sufficiently larger than unity so that the surface wave 
could be excited. 

Fortunately it is possible to engineer such a wave with the proposed structures because 
waves propagating along the sides of the wedge- and pyramidal structures can be considered 
analogous to the surface wave of the chopped medium. A simplified explanation of this 
effect for a pyramidal structure in Fig. HJb) is as follows: At a given frequency the incident 
wave beam excites one of the nanoplates, namely that resonating at this frequency. The 
local field enhancement happens at the edges of the resonant plate and this local field is 
a package of evanescent waves. One of these evanescent spatial harmonics excites at the 
sides of the pyramid a wave propagating from one nanoplate edge to another, i.e., we have 
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FIG. 2. (Color online.) Left - two structures for which the surface wave dispersion curves were 
calculated (chopped interface and solid silver interface). Right - dispersion curves for these struc- 
tures. The dispersion for the chopped interface differs sufficiently from the light line in the whole 
visible range. 



a surface-like wave propagating along the sides of the pyramid to its apex. The field of 
this wave is enhanced compared to the incident wave at the surface of the pyramid and 
decays in the surrounding space the further it is apart from this surface. If the length of the 
nanostructure is properly chosen the local field at the apex of the nanostructure will be be 
enhanced more compared to the sides of the pyramid due to the constructive interference of 
partial surface waves propagating at four sides of the pyramid. The excitation of nanoplates 
which are some distant away from the pyramid base is possible if the whole structure is 
optically rather small and if, for instance, the thickness of silver plates does not exceed 
10 — 15 nm. So thin silver plates are nearly transparent in the visible. 

The operation of the planar structure shown in Fig. EJ^b) differs from the pyramidal 
case only by the dimensionality. The similar edge wave propagating along the sides of the 
tapered planar waveguide can be excited due to the same physical mechanism. In this case 
the requirement of the small optical size of the structure can be strongly softened compared 
to the pyramidal structure. Slender metal nanobars reflect the incident wave beam weakly. 
Therefore the waveguide mode of the input waveguide can be matched with the mode of 
the tapered metamaterial waveguide. We called the eigenwave responsible for the operation 
of both planar and pyramidal structures the edge wave since its electric field is mainly 
concentrated at the nanobar or nanoplate edges. In order to optimize the tapered structure 
the edge mode has been analyzed for non-tapered periodic waveguides. 
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Direction of the periodicity 



FIG. 3. (Color online.) A unit cell of the periodic plasmonic waveguide. 
B. Dispersion of the edge mode 

We numerically studied the dispersion properties of a planar periodic plasmonic waveg- 
uide whose unit cell is shown in Fig. El In this waveguide silver nanobars are sandwiched in 
between dielectric nanoplates. The dispersion characteristics of such a waveguide were stud- 
ied using the HFSS simulator whereas the complex permittivity of the silver corresponded 
to the known experimental data^i and the permittivity of the dielectric was set equal to 
2.2. The dispersion diagrams for different values of the x— dimensions (keeping y— and 
z— dimensions fixed) are shown in Figs. 11(a) and (b). The amplitude of the electric field of 
the different modes at different frequencies are shown in the insets of these Figures. The 
edge modes shown in Fig. 11(a) propagate on the horizontal side of the waveguide whereas 
the modes depicted in Fig. @|b) propagate on the vertical side of the waveguide. 

We see that the modes propagating along the horizontal side of the waveguide are less 
affected by the variation of the x— dimension of the waveguide whereas the dispersion of the 
modes propagating along the vertical side of the waveguide is clearly affected by any change 
in the x— dimension of the waveguide. For our purposes the mode which is less affected by 
the tapering of the waveguide in the x— direction is more favorable. This makes possible 
the tapering of the waveguide and a large operational frequency band. The apex of the 
waveguide can be reduced to sub-wavelength dimensions without sacrificing too much in 
the matching between the wider waveguide sections. The operational band of the wedge is 
determined by the frequency band over which the dispersion curve of the mode on widest 
and narrowest waveguide sections overlap. Clearly, we can see that for the vertical edge 
mode the operational band hardly pushes to the optical regime. For the horizontal edge 
mode on the other hand the operational band covers the whole optical regime. Further, for 
both modes the most part of the electric field energy is concentrated near the surface of 
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FIG. 4. (Color online.) Dispersion diagrams of the modes (a) - on the horizontal edges and 
(b) - on the vertical edges for the plasmonic waveguide with different x— dimensions. The other 
parameters read: y— dimension is 70 nm, z— dimension is 35 nm, and t = 100 nm. The period of 
the waveguide is 135 nm. The amplitude distribution of the electric field over the unit cell at the 
frequency denoted by the black circle is shown in the insets. 

nanobars, outside their bulk. Therefore the energy absorption in these modes is not so high. 
C. Plasmonic resonances of nanoparticles 

One can roughly estimate the resonance frequencies of an individual nanobar (or nanoplate) 
replacing it by a slender (or oblate) ellipsoid. For optically small ellipsoids one can use the 
following condition for the internal electric field E ; :— 



Ei 



-E e -)• oo, 



(1) 



e e + Ni (ei - e, 

where e e is the relative permittivity of the surrounding media, E\ is the relative permittivity 
of the ellipsoid (in our case it is the complex permittivity of silver), and Ni is the depo- 
larization factor along i — x,y,z coordinate axes. For the oblate ellipsoids (spheroids) the 
depolarization factors read:^ 



1 + e 2 

N z = — - — (e — arctan (e)) , 
N x = N y = ±(l-N z ), 



(2) 
(3) 



where eccentricity e = \J a 2 /c 2 — 1. 

The resonance frequency becomes then 

/. = ^ , "* (4) 

In these estimations we neglect losses and the permittivity of silver is assumed to obey the 
Drude dispersion 

= (5) 



where uj p = N q 2 / men is the plasma frequency, N is the charger density, q is the electron 
charge, and m is the electron mass. Naturally, the estimated resonance frequencies calculated 
using the above equations are not quite the ones we have for the rectangular nanoparticles, 
moreover those in presence of dielectric spacers and those coupled with one another through 
the waveguide mode. Nevertheless, this rough approximation gave us a valid estimation of 
the resonance frequency range and of the ratio between different resonance frequencies. 



III. SIMULATED PERFORMANCE OF METAMATERIAL STRUCTURES 



A. Planar nanotip 

We simulated the performance of the planar MMNT in HFSS. A structure comprising 
10 coplanar silver nanobars having thickness d = 70 nm (dimension in y — direction), length 
along the waveguide d = 35 nm (dimension in z — direction) , and length across the waveguide 
(x-dim.) reducing from 256 nm to 64 nm, i.e., following the notations in Fig. [TJ a — 32 nm 
and the length of last 6 nanoplates decreasing with step b = 32 nm. The structure was 
excited with a wave port located at the input. The dielectric spacers had length in z- 
dimension (i.e., p — d) t = 100 nm and same dimensions in x- and y-direction as the silver 
nanobars. 

In Figs. E^a)-(i) we show the distribution of the absolute value of the electric field ampli- 
tude normalized to that calculated in absence of the plasmonic structure. The relative field 
distribution is depicted for the central horizontal plane of the MMNT. 

In all these figures the typical field distribution of the edge wave is clearly seen in the 
narrow part of the MMNT. The field is locally enhanced along the contour of the structure 
until the distances of 30 — 70 nm away from this contour. The amplitude enhancement factor 
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FIG. 5. (Color online.) The normalized amplitude of the electric field at (a) 450 THz, (b) 470 THz, 
(c) 490 THz, (d) 510 THz, (e) 530 THz, (f) 550 THz, (g) 570 THz, (h) 590 THz, and (i) 610 THz. 

is equal 1.5 .. . 7.5 (depending on the frequency). Within the frequency range 450 — 530 
THz the field decreases fast at larger distances from the contour. Within the frequency 
region 530 — 590 THz the field is enhanced not only along the contour. The constructive 
interference of two edge waves observes at these frequencies and leads to a hot spot located 
near the apex. The amplitude enhancement within this hot spot attains 7. . .8 (i.e., the 
intensity enhancement is 50 . . . 60). Naturally, even a so large local field enhancement is not 
comparable with that achievable in self-similar nanolenses.— However, the combination of a 
rather large field enhancement with a relatively large hot spot (its diameter is 30 ... 50 nm), 
a stable location of the hot spot in a very broad frequency band makes our structure very 
attractive, especially if we compare the fabrication costs of our structure to that of a self- 
similar nanolens. At 590 — 610 THz the field enhancement along the contour disappears but 
the hot spot at the apex keeps. At these frequencies the attenuation of the mode becomes 
significant. 
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FIG. 6. (Color online.) Logarithmic distribution of the field intensity over the plane of axial cross 
section of the pyramidal nanostructure, log(I/io). From upper left to down right the frequencies 
read: 862 THz, 680 THz, 590 THz, 526 THz, 462 THz, and 429 THz. 

The main conclusion from Figs. [5] is that the regime which offers the local field enhance- 
ment in front of the end of the MMNT refers to the whole frequency region 450 — 610 THz 
which covers 46 percent of the visible range. In this work we do not consider the matching 
of our nanostructure with a conventional optical waveguide. This problem will be addressed 
in our next paper. 



B. Pyramidal nanotip 

A systematic study of the pyramidal MMNT is not yet finished. In the present paper we 
show results which will be hopefully improved in the next future. Geometrical parameters 
selected below are based only on analytical calculations. Simulations of the nanopyramid 
were done using the FDTD codes.-^^- Two main differences in the operation of the pyramidal 
structure from the planar one are as follows: 1) the multi-frequency operation instead of 
the broadband one since the resonances of the adjacent nanoplates do not overlap, 2) larger 
local enhancement in the hot spot. 
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In simulations whose results are depicted in Fig. [6] the whole space around the silver 
nanoplates was filled with uniform dielectric e = 2.2. Thickness of silver nanoplates in this 
case was 5 nm as well as the thickness of dielectric spacings between the plates. In Fig. [6] 
we show logarithmic field distributions log(///o), where / is the field intensity in presence 
of nanoplates and Io is the field intensity in their absence. The distributions are shown in 
the plane of the pyramid axial cross section. This plane is parallel to the plane in which the 
incident electric field is polarized. The incident wave beam illuminates the pyramid base 
and is directed along its axis. 

The structure resonates at 9 frequencies (385, 429, 462, 526, 588, 632, 681, 769, 862 
THz). Each nanoplate has 2 plasmon resonances except the smallest plate whose shape 
is close to cubic and which has therefore one resonance. For five higher frequencies we 
obtained the hot spot at the structure apex that we treat as the result of the constructive 
interferences of edge waves in this spatial region. For four lower resonance frequencies the 
interference of edge waves near the apex is destructive since the field distribution has a 
minimum there. However at these four frequencies we obtain nearly the same hot spot 
between the second or third nanoplates (counted from the apex) because these frequencies 
correspond to resonances of these two plates in two orthogonal polarizations. So, with this 
structure we obtain one rather stable hot spot for five discrete resonance frequencies and 
another stable hot spot at other four frequencies. We also simulated the structure with 
same 5 nm dielectric spacing between nanoplates but thickness of nanoplates was equal 2 
nm. In such a structure similar situation observes with the location of hot spots but the 
field enhancement in them is significantly larger than that in Fig. [6j 



IV. CONCLUSIONS 

In this paper we suggested, theoretically estimated and simulated using the HFSS package 
and our own FDTD code new tapered nanostructures which is our terminology refer to 
metamaterial nanotips. These nanotips are formed by rectangular silver nanoparticles, either 
nanoplates or nanobars alternating with corresponding dielectric spacers which are arranged 
either in a planar structure or in a pyramidal one. The structures are excited by a wave 
beam incident from the structure bases. Subwavelength spatial regions are formed around 
the structure and/or in front of its apex in which the electric field is locally enhanced. In 
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the space surrounding the structure the field is significantly smaller than the incident one. 
The advantage of our structure compared to known ones is that the subwavelength spatial 
region of the enhanced field is formed nearly the same at different frequencies. For the 
planar structure the frequency stability holds over a very broad frequency range. For the 
pyramidal one this is the stability for several discrete frequencies. We believe that these 
advantageous properties of our structures and their practical feasibility will not only open 
new doors in field-enhanced schemes of nanosensing and microscopy, they will be used in 
prospective nanophotonic applications. 
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